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Long-distance QCD effects in Bd,s — ► £ + £~"y decays are analyzed. We show that the contribution 
of the light vector-meson resonances related to the virtual photon emission from valence quarks of the 
B-meson, which was not considered in previous analyses, strongly influences the dilepton differential 
distrubution. Taking into account photon emission from the 6-quark loop, weak annihilation, and 
Bremsstrahlung from leptons in the final state, we give predictions for dilepton spectrum in Bd, s —> 
£ + £~j decays within the Standard Model. 



PACS numbers: 13.20.He,12.39.Ki,12.39.Pn 
O . I. INTRODUCTION 

o ■ 

Rare radiative leptonic Bd, s — > £ + £ 7 decays are induced by the flavour-changing neutral current transitions 
b — > s, d. In the Standard Model such processes are described by penguin and box diagrams and have small probabilities 
f~\ ■ 10" 8 - 10" 15 (see e.g. Q). Processes with such small branching ratios cannot be observed at the running machines 
such as Tevatron, BaBar and Belle. The decays Bd, s — > and B c i, s — ► fi + fi~"f will be studied at LHC with the 

detectors ATLAS, CMS, and LHCb 0. These decays are perspective candidates for a search for physics beyond the 
\ Standard Model, therefore reliable theoretical predictions for these decays are of great interest. 
The effective Hamiltonian describing the b —> q (q = d,s) weak transition has the form |3| 

^^^xVawo.w, (i.i) 



x 



,— ' ■ Gf being the Fermi constant, Cj the scale-dependent Wilson coefficients, and Oi the basis operators. For B decays 
fi ~ 5 GeV is a convenient choice. The amplitudes of the basis operators between the initial and final states may 
■ be parametrised in terms of the Lorentz-invariant form factors. For radiative leptonic decays several different basis 
i-G ' operators contribute, and respectively, one encounters several different types of form factors. The latter contain 
Oh! nonperturbative QCD contributions, and therefore their calculation poses one of the main problems for theoretical 
O ' analysis of B — > £ + £~^ decays. 

The B — > £ + f~7 decays have been studied in several papers 0, 0, 0, S where transition form factors have 
i-C ■ been analysed and partial widths, photon energy spectra, dilepton mass spectra, and charge asymetries have been 
| calculated. It turns out however that an important contribution related to hadron resonance in the £ + £~ channel 
was not taken properly into account. Our goal is to account for this contribution and to give reliable predictions for 
dilepton mass spectrum and decay rates. 1 

The paper is organized as follows: In Section III Al we discuss the dominant contribution to the decay amplitude 
related to the transition ("f\H c ff(b — > q£ + £~)\B). In Section Til Bl we study the contribution related to the B — > £ + £~ 
transition induced by the photon penguin operator {£ + £~\da^ v q\B). For a proper description of this process it is 
necessary to take into account light meson resonances which emerge in the physical decay region. Making use of the 
vector meson dominance we relate these contributions to the B — > Vj transition form factors. Weak annihilation is 
discussed in Section III CI This process is known to be suppressed by 1/ mj compared to the photon emission from 
the _B-meson loop. Nevertheless, it gives sizeable contribution to the dilepton mass spectrum in the region of small 
dilepton momenta. In Section III Dl we consider Bremsstrahlung from leptons in the final state. Its contribution to 
the cross section is proportional to (m^/Afs) 2 and is thus essentail mainly for r leptons in the final state. In Section 
IIIII we give numerical predictions for the branching ratios and dilepton mass spectra based on our detailed analysis of 
the form factors. 



1 We do not consider any asymmetries since they cannot be studied at LHC because of a small expected number of events 0. 



Fig. 1: Diagrams contributing to Bd — » £ + £ 7 discussed in section III Al Dashed circles denote the b — > d"f operator Oj 1 . Solid 
circles denote the b — > d£ + £~ operators Ogv and Oioav ■ 



II. THE B -> i+r-y AMPLITUDE 

A. Direct emission of the real photon from valence quarks of the B meson 

The amplitude of the process when real photon is directly emitted from the valence b or d quarks, and the £ + £~ 
pair is coupled to the penguin is described by diagrams of Fig. ^ This amplitude is equal to { n /\H^ M 1 \B) with 2 
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The coefficient Cm/{p,,q 2 ) includes long-distance effects related to cc resonances in the g 2 -channel, q the momentum 
of the £ + £~ pair The Ci^ part in Eq. I|2.1|l emerges from the diagrams in Fig. ^with the virtual photon emitted 
from the penguin 
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The B — > 7 transition form factors of the basis operators in (|2.1|) are defined according to 

F A (q 2 ) 



(2.2) 



(l(k, e)\dn/^ 5 b\B{p)) = i e e* (g m pk - p a k,j) 
( 7 (fc, e)\drt»b\B(j3)) = ee*£ K ,?A 
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(2.3) 



(l(k, e)\da^j 5 b\B(p)) (p - k) v = ee* a [g^ a pk - p Q fc M ] F TA (q 2 , 0), 
{'y{k,e)\da llv b\B{p)){p-k) u = i e e* e^^p^kr, F TV (q 2 , 0). 

We treat the penguin form factors Ftv,ta(Qi , ?I) as functions of two variables: q\ is the momentum of the photon 
emitted from the penguin, and qi is the momentum of the photon emitted from the valence quark of the B meson. 
Usually one denotes FTv,TA(q 2 ,0) = Fxv,TA(q 2 )- The form factors FA,v,TA,Tv{q 2 , 0) were studied in detail in 01 . 
The parametrizations for these form factors from £| which satisfy all known constraints in the limit m b — ► 00 0, |n| 
will be used in our analysis. 

The amplitude corresponding to diagrams of Fig 1 takes the form Q| : 
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2 Our notations and conventions are: 7 s = i7°7 1 7 2 7 3 , <r MI , = s[7f»,7»>]) e 0123 = —1, e a bcd = e af3fiv aa b@ c^d" , e = 



Fig. 2: Diagrams contributing to Bj — > £ + £ 7 discussed in section Hi Bl Dashed circles denote the b — > d~f operator O-j-^ 



B. Direct emission of the virtual photon from valence quarks of the B meson 

This process is described by diagrams of Fig. 2: the real photon is emitted from the penguin, whereas one of the 
valence quarks directly emits the virtual photon which then goes into the final pair. The corresponding amplitude 



A^' has the same structure as the CV 7 part of the amplitude A y ^' with Fta,tv(q 2 , 0) replaced by -Fta,tv(0, q 2 )- The 
form factors Fta,tv (0, q 2 ) at the necessary timelike momentum transfers are not known. The difficulty with these 
form factors comes from neutral light vector mesons, p° and u> for B c i decay and (f> for B s decay, which appear in the 
physical B — > jl + £~ decay region. These resonances emerge in the amplitude of the subprocess when the photon is 
emitted from the light valence d or s quark. We obtain the form factors Fta,tv(®> I 2 ) f° r <? 2 > using gauge- invariant 
version of the vector meson dominance This allows us to express the form factors Fta,tv(0, q 2 ) in terms of 
the B — ► V transition form factors at zero momentum transfer and leptonic constants of vector mesons 



F T v,TA{0,q 2 ) = F TVtTA (0,0) - 2 fv9+~ 
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My and Fy being the mass and width of the vector meson resonance. The B 
according to relations 



V transition form factors are defined 



(V(q,e)\da^b\B(p)) = ie* a e^ [g^ v (k 2 )g a p(p + q)i + g 
The leptonic decay constant of a vector meson is given by 

(0\d^d\V{e,p))=s fi M v f v 



B^V 



(k*)g a pki + g^ v (e)p a pPq 



(2.6) 



(2.7) 



We shall use the form factors g B ~* v calculated in using the relativistic dispersion approach [Tfij . 

Notice that the resonance contribution discussed in this Section is similar to the resonance contribution in the 
Wilson coefficient Cg V which emerge due to 4-quark operators. In the case of Cgy vector mesons containing cc and 
uu pairs contribute. 



Weak annihilation 



The weak annihilation amplitude A WA is given by triangle diagrams of Fig 3. One should take into account u and 
c quarks in the loop. The vertex describing the bd — > UU transition (U = u, c) reads 



H, 



B^UU 



off 



Gf 

71 



01 VuhVZ A d 7/1 (l - i 5 )b U-y^l - -ftp, 



(2.8) 




Fig. 3: Weak annihilation diagrams discussed in section 111 CI Triangle diagrams with u and c quarks in the loop should be 
taken into account. 
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with a\ = Ci + C2/N c , N c number of colors |16|, For N c = 3 one finds ai = —0.13. We now have to take 



The C/C/ contribution to this amplitude can be written as 

G 



B^UU 



B). 



G 77 (M 2 B ,k 2 = 0,q 2 \m 2 u ) 



V2 " q 

where the form factor G(p 2 , k 2 , g^m 2 /) is defined as follows 01 

(YikWiQWiU-y^Um = e 2 e* a (k)e^(q)e a p kq G 11 (k 2 ,q 2 ,p 2 \m 2 u 



(2.9) 



(2.10) 



(2.11) 



a i — ■ For c- 



For massless u-quark in the loop, axial anomaly ^| fixes the form factor G 77 (p 2 , k 2 , q 2 \Q) = ~ 2Nc ^"^ ■ 
quark there is additional g 2 -dependent contribution given by the amplitude m c (7*(fc)7*(<7)|c75c)|0) ~ m 2 c /Mg, which 
contains ip and ip' resonances at q 2 > 0. The latter contribution is numerically negligible compared to contributions 
discussed in the previous sections for all q 2 in the reaction of interest. Therefore, we have 



.4 



WA 



^|a cm e ai{V ub V* d + V cb V* d }^e^ qk ^ l+^^t 



(2.12) 



The anomalous contribution is enhanced at small q 2 , but even here it is suppressed by a power of a heavy quark mass 
compared to the contributions discussed in the previous sections [ill Il9| . 



D. Bremsstrahlung 

Fig. 4 gives diagrams describing Bremsstrahlung. The corresponding contribution to the B — ► £ + £~~/ amplitude 
reads 



i Brcms 



Gf «em T ,*, 



h 



' B 7f ~2T V * dVtb Mb ^ CwA ^ lijp2) 



( 7 e*) ( 7 p) (75*) 



i — rhj 



75 £{-pi)- 



(2.13) 





Fig. 4: Diagrams describing photon Bremsstrahlung. Solid circles denote the operator Oioa- 
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III. THE B -» i + £-j CROSS-SECTION 



The amplitudes discussed in Sections III AIII (Jl have the same Lorentz structure, whereas the structure of the 
Bremsstrahlung amplitude ( Section III D|) is different. Therefore, it is convenient to write the cross-section as the sum 
of three contributions: square of the amplitude A 1+2+WA which we denote square of the amplitude ^4 Brcms which 
we denote T( 2 \ and their mixing, denoted as I^ 12 ); 



G 



2 ^3 
F em 



Ml 



\Vt b v; q \ 



x 2 Bo (s, t) + x £ (s, t) B x (s, t) + C 2 (s, i) B 2 (s, t) 



d 2 T^ _ 

dsdt ~ 2 W 7T i 

Bo (J, t) = (S + Am 2 ) (Fx (a) + F 2 (a)) - 8mj |C 10 a(m)| 2 (K (? 2 ) + F 2 A (q 2 )) , 
B^sJ) = 8 [sF v {q 2 )F A {q 2 )Re(cl{/*{^q 2 )C W A^)) 

+ m b F v {q 2 ) Re (C 7 %( M ) F^ A (q 2 ) C WA (jj)) + rh b F A {q 2 ) Re (C^ifi) F^ v (q 2 ) C WA (fi))] 
B 2 (a, t) = a (Fx (8) + F 2 (J)) , 



(3.1) 



M*> = (|C9v(M,3 2 )| 2 + \C 10 a(v)\ 2 ) F 2 (q 2 ) + (-^) \C 7l ^)F TV {q 2 ) 



2m b 



+ ^ F v (q 2 ) Re (c 7 » F TV (q 2 ) C e 9 {/*(», q 2 )) 



Fr, 



(s) = (|C 9 W,M)| 2 + \C WA {n.)\ 2 ) F 2 A {q 2 ) + (^) \C 7 ^)F TA {q 2 ) 
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+ ^ F A (q 2 ) Re ( CV 7 (/i) F TA (q 2 ) C 9 e p(/x, 9*) 
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ds dt 

rf 2 r< 12 ) _ _^ 

dsdt ~ 2 10 7T 4 

2 x m(, 
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M B 
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s + x 2 /2 



(u - m 2 e )(t - mf) \(u - mj) (t - mf) 



x mi 



tg\ M 



B q 4 2 



(it - m 2 )(t - mf) 



(3.2) 
(3.3) 



Re(C* WA {n)C 7l (p)F TV (q 2 ,Q)) + x F v (q 2 ) Re (C* 10A (^C^(^ q 2 )) + £(a,t)F A (q 2 ) \C WA (n)\ 



Here 



(p - ky 

M 2 B 



t = 



(g - Pi) 
Ml 



with s + t + u = 1 + 2m 2 , m 



:j - mf/M§, m b = m b /M B and 



(p - ggT 

M 2 



£ (s, t) — u — i. 



x \/l — 4m 2 /s ' 

In the above formulas the complex form factors Ftv,ta and defines as follows 

16 V ub V* d + V cb V* d ax / fl 



(3.4) 



(3.5) 



^rvfaT) - F Ty (^,0)+ J F T v(0,^)- 

^(g 2 ) = F TA {q 2 ,0) + F TA {0,q 2 



VtbV* d C 7l m b ' 



(3.6) 



The expressions (I3.2|) and (|3.3|l contain the infrared pole which requires a cut in the energy of the emitted photon. 
Clearly, the contribution of the pole is proportional to the lepton mass. 
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IV. NUMERICAL ANALYSIS 
A. Parameters 



For numerical estimates we use the following values: 

• The Wilson coefficients at n = 5 GeV corresponding to C%(Mw) = — 1: 

d{5GeV) = 0.241, C 2 (5 GeV) = -1.1, C 7l (5GeV) = 0.312, C 9V (5GeV) = -4.21, C 10A (5GeV) =4.64 

• The CKM matrix elements \V^V td \ = (8.3 ± 1.6)10- 3 , \V t * b V ts \ = (4.7 ± 0.8)10" 2 H3- 

• The B-meson lifetimes r(B d ) = 1.536 ± 0.014 ps and t(B s ) = 1.461 ± 0.057 ps HJ. 

• For the Bd — > 7 form factors we use parametrizations from £| 



F l (q 2 ) = ft 



!bM l 



E 1 = 



Mr 



'Ai + W 2 
with i = V, A, TV, TA and the parameters (3 and A given in Table 



1 - 



(4.1) 



Table 1: Parameters for the form factors in Eq. <|4.1jl from 



Parameter 


Fv 


Ftv 


Fa 


Fta 


(5 (GeV" 1 ) 
A (GeV) 


0.28 
0.04 


0.30 
0.04 


0.26 
0.30 


0.33 
0.30 



• The form factors g+~* v at q 2 = for — > p° , — > w transitions and £? s 
-B s — > 7 and the i?d — > 7 form factors are taken to be equal to each other. 



are given in Table [3 The 



Table 2: The leptonic decay constants fv and the B — > V transition form factors g+^ v at q 2 = 0. 





P° 




<$> 


fv (MeV) 


154 


45.3 


-58.8 


fff^(O) 


0.19 


-0.19 


-0.38 



B. Results 

The calculated dilepton differential distributions are shown in Fig. EJfor different leptons in the final state and for 
different values of the photon energy cut - the minimal photon energy in the i?-meson rest frame -E^in- Table [31 gives 
the dependence of the integrated Bd, s — * l + £~j decay rate on E^ in . A particular choice of £^ in depends on the 
energy resolution in a specific experiment: namely, E^ in — 80 MeV and E^ in — 20 MeV in Table 01 correspond to 
the expected accuracy of the B-meson reconstruction of ATLAS and LHCb, respectively. 



Table 3: The Bd,s — > £ + £ 7 decay rates as functions of the minimal photon energy The region of the J/tp and ip' 

resonances 0.33 < s < 0.55 was excluded from the s integration, that corresponds to the experimental procedure adopted at 
LHC HEl 



rat 


m e 




m T 


El, in (MeV) 


20 


50 


80 


20 


50 


80 


20 


50 


80 


Br (B d -» e+Tj) x 10 10 [This work] 
Br (B s -> £ + r~f) x 10 9 [This work] 


3.95 
24.6 


3.95 
24.6 


3.95 
24.6 


1.34 
18,9 


1.32 
18.8 


1.31 
18.8 


3,39 
11.6 


2.37 
8.10 


1.87 
6.42 


Br (B d -> £+£-7) x 10 10 [5J 
Br (B„ -» i+e-j) x 10 9 [5] 


1.01 
3.30 


1.01 
3.29 


1.01 
3.29 


0.66 
2.16 


0.62 
2.06 


0.61 
2.00 


3.39 
11.6 


2.38 
8.15 


1.88 
6.47 


Br (B s -> £ + r-y) x 10 9 [6j 


20 


20 


20 


12 


12 


12 









Fig. 5: Dilepton g 2 -spectrum in B s — » e + e 7 (left), B s —* fi + /x 7 (central), and B s — > r + r 7 (right) for different photon 
energy cuts: E^ in = 20 MeV (solid), E^ in = 50 MeV (dashed), EP min = 80 MeV (dotted). 



Our results for the integrated decay rate in Table 03 are close to the results reported in Q . Notice however that 
in [|| resonance contributions were not considered and the B —> 7 form factors corresponding to the limit nib — * 00 
were used. As shown in Q, corrections to the asymptotic formulas in B — > 7 form factors at large q 2 are rather 
large, therefore the differential distributions calculated here differ strongly from those of 6]. 3 Our results both for 
the integrated and the differential distributions in the B — > e + e~7 and B — > /i + /i~7 decays in the region q 2 < 2 GeV 
differ strongly from the results of |5J where the contribution of the direct virtual photon emission from valence quarks 
of the B-meson (Section fll B(l was not taken into account (Figs. El Cl- 




io 10 2 10 1 0.5 0.6 0.7 0.B 0.9 

q ! /M', q*/UC 



Fig. 6: Dilepton q 2 -spectrum in Bd decays: Bd — > e + e~7 (left), Bd — * fJ-~lu-~J (central), and Bd — > r + r~7 (right) for 
= 20 MeV. Solid line - our result; dashed line - result corresponding to .5|, where the contribution of the direct virtual 
photon emission from valence quarks of the B meson (Section |TT3 was not taken into account. 



V. CONCLUSIONS 

We studied the Bd, s — * £ + £~"f decays taking into account photon emission from the 6-quark loop, weak annihilation, 
and Bremsstrahlung from leptons in the final state. Special emphasis was laid on long-distance QCD effects in 
B — ► £ + £~7 decays related to the photon emission from the 6-quark loop. The contribution of light vector-meson 



3 Let us point out that in distinction to the J /ip and ip 1 , the region of light vector resonances will not be excluded from the experimental 
analysis. That is because for small values of s photons have sufficient energies to be registered in the electromagnetic calorimeter. 
Moreover, the region of small s gives the main contribution to the B — » 7 signal at LHC. For B s — ► e+e _ 7 and B 3 — * fi+fi~j 

decays, the 0-meson resonance contribution leads to a strong enhancement of the full spectrum compared to the non-resonance one. 
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10 J 0.5 0.6 0.7 0.8 0.9 

q 2 /M* B q'/W, q ! /Mj 



Fig. 7: Dilepton q 2 -spectrum in B a decays: B a — > e + e 7 (left), B s — > fi 7 (central), and B s — > r + r 7 (right) for = 20 

MeV. Solid line - our result, dashed line - [f|, dotted line - [g. 



resonances related to the direct virtual photon emission from valence quarks of the B-meson, not taken into account 
in previous analyses, was shown to be essential for a proper description of the process. In particular, this contribution 
affects stongly dilepton mass spectra in processes with light leptons in the final state. Reliable predictions for dilepton 
mass spectrum in Bd. s ~> £ + £~j decays within the Standard Model were given. 
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